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SECTION  I 


INTRODUCTION 


Cables  containing  many  wires  have  been  widely  used 
for  signal  transmission  within  large  systems.  Recently, 
the  effects  of  electromagnetic  pulse  (EMP)  have  drawn  much 
attention  to  the  studies  of  effective  shielding  techniques 
for  such  multiconductor  cables.  A practical  shielding 
structure  for  these  cables  is  a grounded  conducting  tube,  as 
shov/n  in  Tigure  1.  In  order  to  study  the  mutual  coupling  of 
the  wires  in  the  presence  of  the  grounded  shield,  the  per- 
unit-length  capacitance  matrix,  which  determines  the  "quasi-TEM" 
transmission  line  solutions  (ref.  1,  2,  3),  must  be  obtained. 

The  purpose  of  this  report  is  to  demonstrate  a numerical 
method  which  calculates  the  per-unit-lencth  capacitance  matrix 


1.  F.Y.  Chang,  "Transient  Analysis  of  Lossless  Coupled  Trans- 
mission Lines  in  a Nonhomogeneous  Dielectric  Medium," 

IEEE  Trans,  on  Microwave  Theory  and  Techniques,  Vol.  MTT-18, 
No.  9,  pp.  616-626,  September  1970. 

2.  C.R.  Paul,  "On  Uniform  Multimode  Transmission  Lines," 

IEEE  Trans,  on  Microwave  Theory  and  Techniques,  Vol.  MTT-21, 
No.  8,  pp.  556-558,  August  1973. 

3.  C.R.  Paul,  "Useful  Matrix  Chain  Parameter  Identities  for 
the  Analysis  of  Multiconductor  Transmission  Lines," 

IEEE  Trans,  on  Microwave  Theory  and  Techniques,  Vol. 

No.  9,  pp.  756-760,  September  1975. 


MTT-23 , 


(ref.  4,  5,  6)  of  the  dielectric-coated  wires  in  Figure  1.  The 

i^^  wire  has  a conducting  core  of  radius  a^  , the  outer  radius 

of  the  dielectric  layer  is  , the  relative  dielectric 

constant  is  , with  the  value  of  i=l,  2,  ...,n.  The 

total  number  of  wires  is  n . The  center  of  the  i^^  wire 

is  • The  shielding  tube  has  its  center  located 

at  the  origin  (x„  = 0 , y =0)  with  a radius  of  c . 

c c 

Theoretically,  the  problem  can  be  solved  by  the  con- 
formal mapping  technique  (ref.  4,  5).  However,  the  dielectric 
layer  surrounding  the  wires  makes  the  mapping  very  complicated. 
It  is  also  possible  to  solve  the  problem  using  the  method  of 
successive  images  (ref.  4,  7).  Again  the  process  of  images 
involving  the  dielectric  coatings  is  not  simple,  and  the 
imaging  is  very  tedious  if  many  wires  are  involved. 


4.  R.  Plonsey  and  R.  Collin,  Principles  and  Applications  of 
Electromagnetic  Fields,  McGraw-Hill  Co. , New  York  (1951) , 
p.  71. 

5.  E.  Weber,  Electromagnetic  Fields  - Theory  and  Applications, 
Vol.  1 - Mapping  of  Fields,  John  Wiley  and  Sons,  New  York 
(1950) , pp.  139,  223,  226,  379. 

6.  J.C.  Clements,  C.R.  Paul,  A.T.  Adams,  "Computation  of  the 
Capacitance  Matrix  for  Systems  of  Dielectric-Coated 
Conductors,"  IEEE  Trans,  on  Electromagnetic  Compatibility, 
Vol.  EMC- 17,  No.  4,  pp.  238-248,  November  1975. 

7.  M.P.  Sarma  and  W.  Janischewskj , "Electrostatic  Field  of 

a System  of  Parallel  Cylindrical  Conductors,"  IEEE  Trans. 
on  Power  Apparatus  and  Systems,  Vol.  PAS-88,  No.  7, 
pp.  1069-1079,  July  1969. 
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Recently,  the  integral  equations  (ref.  8)  and  the  method 
of  moments  (ref.  6,  9,  10,  11,  12)  have  been  applied  to  solve  < 

the  capacitance  matrix  for  a system  of  multiconductor  wires. 

The  accuracy  and  numerical  stability  of  the  point-matching  tech- 
nique employed  in  references  (6)  and  (12)  are  sensitive  to  the 
points  chosen  for  matching  the  boundary  conditions.  In  addition, 
there  is  no  systematic  way  of  selecting  the  match  points  when 
the  wires  are  randomly  scattered. 

In  this  report,  we  shall  provide  a theoretical  back- 
ground for  the  "multi-series  expansion"  of  the  potential  exterior 
to  the  dielectric  coatings.  The  coefficients  in  this  expansion 
will  be  solved  by  the  "weighted  least-square  method"  and  the 
per-unit-length  capacitance  matrix  is  then  obtained  directly 
from  the  expansion  coefficients. 

3.  P.C.  Chestnut,  "On  Determining  the  Caoacitances  of  Shielded 
Multiconductor  Transmission  Lines,"  IEEE  Trans,  on  Microwave 
Theory  and  Techniques,  Vol.  MTT-17,  No.  10,  pp.  734-745, 

October  1969. 

9.  R.F.  Harrington,  Field  Computation  by  Moment  Methods, 

MacMillan  Co.,  New  York  (1968). 

10.  A.T.  Adams,  Electromagnetics  for  Engineers,  Ronald  Press 
Co.,  New  York  (1972). 

11.  M.P.  Sarma,  "Application  of  Moment  Methods  to  the  Computation 
of  Electrostatic  Fields  - Part  I - Parallel  Cylindrical 
Conductor  Systems,"  IEEE  PES  Summer  Meeting,  San  Francisco, 

California,  July  9-14,  1972. 

12.  A.E.  Feather  and  C.R.  Paul,  Computation  of  the  Capacitance 
Matrices  for  Ribbon  Cables,  RADC-TR-76-101 , '’ol.  II, 

Rome  Air  Development  Center,  NY,  April  1976. 
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SECTION  II 


THEORY 

In  this  section  we  shall  describe  the  field  expansions 
in  all  regions  of  the  multiconductor  cable.  The  "multi-series 
expansion"  will  be  derived  to  represent  the  potential  exterior 
to  the  dielectric  coatings.  The  per-unit-length  capacitance 
matrix  will  be  directly  related  to  the  expansion  coefficients. 

1.  Single  Wire  Solutions 

Consider  the  conductor  i,  as  in  Figure  2a.  The  region 

Ri  (a^  < r^  < b^)  is  filled  with  dielectric  material  with 

dielectric  constant  . Outside  the  coating  is  the  region 

R (r.  > b)  with  dielectric  constant  . According  to  refer 
oil  o 

ence  (4) , a complete  solution  of  the  electrostatic  ootential 

. . in  reaion  R.  is 
1 1 


m=l 


^ '"io  iT^  ^ ^io 


ai  < Ti  < b.  (1) 


In  Figure  2b,  the  solution  of  the  region  R|  , which  has 
dielectric  constant  , is 


m=l 


< b^  (2) 


7 


Hence  the  solution  of  region  , due  to  the  i’'“  con- 

ductor, in  both  Figure  2a  and  Figure  2b  is  the  same  if  the 
equivalent  charge  density 


^si  3n  t^=b^ 


The  general  expression  of  4;^^  is 


)i  = L ""^i  m$  . ) r 


^ ^iO  ET  ^ Ho 


where  and  ^ are  the  potential  expansion  coefficients 

for  the  exterior  region,  R . The  eauivalence  of  ip  . in 

o ■ oi 

both  Figures  2a  and  2b  induced  by  the  surface  charae  density 
Ogi  is  due  to  the  induction  theorem  as  described  in  ref.  (13) . 

2.  The  Multi-series  Expansions 

If  the  induction  theorem  discussed  above  is  applied  for 

the  conductors  i,  i = 1,  2,  n,  the  solution  exterior 

to  the  dielectric  coatings  of  all  the  wires  of  Figure  1 is 

equivalent  to  that  of  Figure  3.  The  equivalent  surface  charge 

densities  t)er-u.nit-lenqth  are  o o o on  each  di- 

’ si  s2  sn 

electric  surface  and  o „ on  the  surface  of  the  grounded  tube. 

su 

Hence,  by  using  the  superposition  theorem  (ref.  4),  the  complete 
expansion  of  the  potential  is 

13.  R.F.  Harrington,  Time-Harmonic  Electromagnetic  Fields, 
McGraw-Hill,  New  York  (1961). 
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n p 00 


i=l  _ m=l 


-m 


r . 

* "io  *"  57 


m=l 


^ ^oO  H ^ ^oO 


Ti  > bi  , r < c (5) 


where  the  constants  are  all  absorbed  in  ; and 


oO 


^ i = 1 , 2 , . . . , n ; and  r < c . 

In  equation  (5)  the  potential  due  to  the  first  summation  is 

associated  with  all  the  interior  conductors,  whereas  that 
due  to  the  second  sum  is  associated  with  the  shieldinq  tube. 

In  addition,  the  potential  in  the  dielectric  coating 
region  of  conductor  i is  given  in  equation  (1)  as 

30 

Z/  ^ \ / in  . “HI, 

(aimcos  mcfc.  + S.^sin  m<i..)  (r-  + y.^  r^  ) 

m=l 


1 j 

+ a . „ in  — + o . - 
"lO  lO 


where  "i  ^i  ' i = i»  2,  ...,  n . 


[6) 


3.  Boundary  Conditions 

On  each  conducting  surface,  r^  = , the  potential 

is  Vj  . That  is 


11 


[ 


1 D 


in=l 


= V, 


for  0 < ({)j  < 2it 


(7) 


Since  equation  (7)  must  hold  for  all  0 _ (J)  j - 2tt  , we  have 


2m 

"jm  = -^j 


m=  1/  2,  3, 


and  j = 1 , 2 1 


(8) 


and 


3jo  = Vj  ; j - 1,  2,  . . . , n 


(9! 


On  each  dielectric  surface,  potential 

is  continuous  and  the  free  charge  is  zero.  Hence, 


t'i'-i  - ’i'lr.=b.  ■ ° 


(10) 


: J 


and 


L":  3r. 


= 0 


(11) 


3 "i 


where  'ii  and  are  given  in  equations  (5)  and  (6) 


respectively . 

On  the  grounded  tube,  the  potential  is  grounded  to 

zero , 


12 


[■Hr=c  = 0 


(12) 


One  further  condition  is  required  to  furnish  the  solution 
of  the  coefficients.  This  condition  is  obtained  from  the 
conservation  of  electric  charges.  That  is 


2Tr  r 1 ^ r H ^ 

I ••  [S]..,  •■II  ‘‘  4 

n 1=1  o'-  J 


ri=ai 


(13) 


4.  The  Per-Unit-Length  Capacitance  Matrix 

Let  the  total  free  charge  per-unit-length  on  the 


conductor  i be  Q|  . Then 


2Tr 


q:  - - 


/ 1 ar. 

o ^ Jr . =a . 


d<J)^ 


1 ^i 


■2^  ^i  “io  ^i  iT  = 


-2Tr  f-  . a . „ 
1 lO 


(14) 


The  per-unit-length  capacitance  matrix  is  defined  as 


q; 


q; 


q; 


c'  c 


11  "'12' 


C' 


In 


C’  Cl 


21  ''22' 


Cl 


2n 


C'  C'  , 
nl  n2 


C' 


nn 


V, 


V, 


V 


(15) 


13 


J 


Taking 


V. 

1 


if  i 


we  have 


1 volt 
0 


If  1 = : 
otherwise 


(16 


The  complete  per-unit-length  capacitance  matrix  is  obtained 
by  taking  j = 1,  2,  n. 


SECTION  III 


CALCULATION  OF  THE  EXPANSION  COEFFICIENTS 

Plug  and  ' which  are  given  in  equations 

(8)  and  (9) , into  equation  (6) , and  truncate  the  expansion 
series  into  M.  finite  terms.  We  then  have 

/ a^"'  \ 

m=l  ' ^ 


+ cx^Q  ii-H  , i — Ij  2,  •..>  n (171 

Similarly,  summing  the  expansion  series  to  terms  in 

equation  (5) 

n ^i 

= Z Z (n.j„cosm4)i  + 5in,sin  mg.)  r."* 
i=l  m=l 


M, 


r . ] 

+ n £n  — + ) (n  cos  mg  + £ sin  mO)  r 

^iO  b . cm  om 

^ m=l 


.iti 


+ n « S rt 

oO  c oO 


(18) 


The  coefficients  , 3^^^  r i ' -^om 

in  equations  (171  and  (18)  are  unknown  constants  which  remain 
to  be  solved.  The  total  number  of  unknowns  N is 


15 


I 


N = 4 


n 


Mi  + 2Mq  + 2n  + 2 


(19) 


These  unknown  coefficients  can  be  solved  by  using  the  boundary 
conditions  given  in  equations  (10),  (11),  (12),  and  (13). 

It  is  possible  to  use  the  point  matching  technique 
(ref.  12)  to  set  up  the  linear  equations  for  these  coef- 
ficients. However,  the  stability  of  the  linear  equations 
is  very  sensitive  to  the  matching  points  selected.  For 
general  problems,  it  is  necessary  to  employ  a systematic 
way  of  establishing  a stable  system  of  linear  equations. 

From  equation  (10)  , multiply  both  sides  by  l,cos  m(|)j  , 
sin  m(t)^,  m = 1,  2,  ...»  successively,  and  integrate  from 

= 0 to  2t  , we  have 


^TT 

0 

2tt 

cos  m({).  di  . = 0 


(20a) 


(20b) 


j j 


^ dOj  = 0 


where  m = 1,  2,  ...,  M-  and  j = 1,  2, 
n ^ 


lishes  2 ^ M.  + n 


equations . 


, n . This  estab- 


16 


Similarly,  the  condition  given  in  equation  (11) 

n 

suffices  for  another  2 ^ M-  + n linear  equations.  Bv 

j=l  J 

multiplying  equation  (12)  by  1,  cos  m<J)  , sin  mij)  , m = 1,  2,  . 
Mq  successively,  we  have  2M^  + 1 more  equations.  Adding 
the  condition  of  equation  (13) , we  thus  have  a total  of 


N . 4 ^ M. 
i=l 


+ 2Mq  + 2n  + 2 


equations.  The  coefficients  are  then  solved  by  the  Gaussian 
elimination  algorithms  of  the  system  of  linear  equations. 
Once  the  coefficients  are  solved,  the  per-unit-length 
capacitance  matrix  is  obtained  via  Section  II-4. 
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SECTION  IV 
NUMERICAL  RESULTS 

A FORTRAN  program  has  been  written  to  implement  the 
techniques  described  in  this  report.  The  program  is  verified 
by  the  results  of  the  three  cases  given  in  Figures  4a,  4b, 
and  4c.  The  first  case  to  check  the  program  is  the  coaxial 
cable  with  a layer  of  dielectric  constant  t . The  exact 
per-unit-length  capacitance  is 


The  numerical  results  of  all  the  values  of  a,  b,  c,  and  e 
used  in  the  computations  all  agree  with  the  exact  per-unit- 
length  capacitance  to  five-digit  accuracies. 

The  per-unit-length  capacitance  matrix  of  the  two 
conductors  of  Figures  4b  and  4c  are  also  compared  aaainst  the 
exact  per-unit-length  capacitances  (ref.  14).  Five-dicjit 
accuracies  are  obtained  for  all  the  cases  tested.  The  compu- 
tation time  on  a CDC  7600  computer  for  the  two  conductors 
is  about  0.4  second. 


14.  ITT  Reference  Data  for  Radio  Engineers,  Chapter  24, 
Howard  W.  Sams  & Co.,  Inc.,  Indianapolis  (1975). 
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A Coaxial  Cable 


Two  Conducting 
Wires 


(b 


Two  Conducting 
Wires 


o 


Figure  4 


The  characteristic  numbers  of  Figure  5 are  given 

below 

a^  = 1 unit  length  , i=l,  2,  7. 

b^  = 2 unit  lengths  , i=l,  2, 
c = 10  unit  lengths 


<^cl  ' 

^cl^ 

— 

(0,0) 

(^c2  ' 

yc2^ 

= 

(5,0) 

(^c3  ' 

^c3) 

= 

(2. 5, 4. 3) 

^c4’ 

= 

(-2. 5,4.3) 

^^c5  ' 

^cS^ 

= 

(-5,0) 

^c6) 

= 

(-2. 5,-4. 3) 

^^c7  ' 

^c7^ 

= 

(2. 5, -4. 3) 

M^=6  , i=0,  1,  2,  7. 

For  these  values  and  a dielectric  constant  e . = 2 e 
1 = 1,  2,  ...  7,  the  resulting  per-unit-length  capacitance 
matrix  (in  picofarad/meter)  is 


55.40 

-9.106 

-9.106 

-9.106 

-9.106 

-9.106 

-9.106 

-9.106 

48.35 

-10.75 

-0.2226  -0.0336 

-0.2226 

-10.75 

-9.106 

-10.75 

48.35 

-10.75 

-0.2226 

-0.0336 

-0 . 2226 

-9.106 

-0.2226 

-10.75 

48.35 

-10.75 

-0.2226 

-0.0336 

-9.106 

-0.0336 

-0.2226 

-10.75 

48.35 

-10.75 

-0.2226 

-9.106 

-0.2226 

-0.0336 

-0.2226 -10.75 

48.35 

-10.75 

-9.106 

-10.75 

-0.2226 

-0.0336 

-0.2226 

-10.75 

48.35 
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Similarly,  for  wires,  the  per-unit 

length  matrix  (in  picofarad/meter)  is  found  to  be 


74.55 

-12.35 

-12.35 

-12.35 

-12.35 

-12.35 

-12.35 

-12.35 

61.47 

-14.83 

-0.168 

-0.0212 

-0.168 

-14.83 

-12.35 

-14.83 

61.47 

-14.83 

-0.168 

-0.0212 

-0.168 

-12.35 

-0.168 

-14.83 

61.47 

-14.83 

-0.168 

-0.0212 

-12.35 

-0.0212 

-0.168 

-14.83 

61.47 

-14.83 

-0.168 

-12.35 

-0.168 

-0.0212 

-0.168 

-14.83 

61.47 

-14.83 

-12 .35 

-14.83 

-0.168 

-0.0212 

-0.168 

-14.83 

61.47 

The  computation  time  for  this  7-wire  problem  is  about 
7 seconds,  which  is  reasonably  fast  for  most  practical  compu- 
tations . 

The  inductance  per-unit-length  matrix  L'  can  be  readily 
obtained.  Elements  of  L'  are  independent  of  the  dielectric 
materials , and  hence  for  the  configuration  of  Figure  5 with 
e^=£Q,  i=l,  ...,7,  we  have 


c-  l;-=  ^ u 


(21) 


where  is  the  per-unit-length  capacitance  matrix  of  this 

configuration  computed  by  using  the  computer  code.  U is 
the  7x7  unity  matri.x  and  v is  the  speed  of  light  inside  the 
tube.  For  the  geometry  previously  used,  the  per-unit-length 
capacitance  matri.x  (in  picof ara/meter ) is 
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37.0 

-5.96 

-5.96 

-5.96 

-5.96 

-5.96 

-5.96 

-5.96 

34.7 

-6.76 

-0.258 

-0.0516 

-0.258 

-6.76 

-5.96 

-6.76 

34.7 

-6.76 

-0.258 

-0.0516 

-0.258 

-5.96 

-0.258 

-6.76 

34 . 7 

-6.76 

-0.258 

-0.0516 

-5.96 

-0.0516 

-0.258 

-6.76 

34.7 

-6.76 

-0.258 

-5.96 

-0.258 

-0.0516 

-0.258 

-6.76 

34.7 

-6.76 

-5.96 

-6.76 

-0.258 

-0.0516 

-0.258 

-6.76 

34.7 

The  per-unit-length  inductance  miitrix  (microhenrv/meter ) , 
using  equation  (21)  is 


304 

105 

105 

105 

105 

105 

105 

105 

293 

97 

53 

44 

53 

97 

105 

97 

293 

97 

53 

44 

53 

105 

53 

97 

293 

97 

53 

44 

105 

44 

53 

97 

293 

97 

53 

105 

53 

44 

53 

97 

293 

97 

105 

97 

53 

44 

53 

97 

293 

This  method  of  evaluating  L'  enables  one  to  compute 
the  propagation  matrix  y of  a multiconductor  cable  with 
dielectric  coating,  as 


=2 

Y 


C'L* 
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for  the  current  mode  (ref.  15) . 


15.  F.M.  Tesche  and  T.K.  Liu,  Selected  Topics  in  Transmission 
Line  Theo  ry , AFWL-TR-7  7-7  3']  Air  Force' "'Jeapons  Laboratory, 
Kirtland  Air  Force  Base,  NM,  March  1977. 
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SECTION  V 


CONCLUSIONS 

The  per-unit-length  capacitance  matrix  of  a multi- 
conductor transmission  line  enclosed  in  a shielding  tube 
is  evaluated  by  the  method  of  multi-series  expansion. 

Numerical  values  of  typical  examples  are  given. 

The  per-unit-length  inductance  matrix  can  be  obtained  by 
finding  the  inverse  of  the  corresponding  per-unit-length 
capacitance  matrix  with  all  dielectric  coatings  of  the 
conductors  removed.  The  propagation  matrix  of  the  configura- 
tion is  then  given  by  this  inductance  matrix  and  the  per- 
unit-length  capacitance  matrix  computed  for  the  system  with 
dielectric  coatings. 
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